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ABSTRACT 

Context. KIC 10526294 was recently discovered to be a very slowly rotating and slowly pulsating late B-type star. Its 19 consecutive 
dipole gravity modes constitute a series with almost constant period spacing. This unique collection of identified modes probes the 
near-core environment of this star and holds the potential to reveal the size and structure of the overshooting zone above the convective 
core, as well as the mixing properties of the star. 

Aims. We revisit the asteroseismic modelling of this star with specific emphasis on the properties of the core overshooting, while 
considering additional diffusive mixing throughout the radiative envelope of the star. 

Methods. We pursued forward seismic modelling based on adiabatic eigenfrequencies of equilibrium models for eight extensive 
evolutionary grids tuned to KIC 10526294 by varying the initial mass, metallicity, chemical mixture, and the extent of the overshooting 
layer on top of the convective core. We examined models for both OP and OPAL opacities and tested the occurrence of extra diffusive 
mixing throughout the radiative interior. 

Results. We find a tight mass-metallicity relation within the ranges M G [3.13,3.25] M© and Z G [0.014,0.028]. We deduce that an 
exponentially decaying diffusive core overshooting prescription describes the seismic data better than a step function formulation and 
derive a value of /ov between 0.017 and 0.018. Moreover, the inclusion of extra diffusive mixing with a value of log between 1.75 
and 2.00 dex (with in cm^ s"^) improves the goodness-of-fit based on the observed and modelled frequencies by a factor ~ 11 
compared to the case where no extra mixing is considered, irrespective of the (M, Z) combination within the allowed seismic range. 
Conclusions. We find that the inclusion of diffusive mixing in addition to core overshooting is essential to explain the structure in 
the observed period spacing pattern of this star. Moreover, for the input physics and chemical mixtures we investigated, we deduce 
that an exponentially decaying prescription for the core overshooting is to be preferred over a step function, regardless of the adopted 
mixture or choice of opacity tables. Our best models for KIC 10526294 approach the seismic data to a level that they can serve future 
inversion of its stellar structure. 

Key words. Asteroseismology - Stars: oscillations (including pulsations) - Stars: interiors - Stars: evolution - Stars: individual: 
KIC 10526294 


1. Introduction 

Intermediate and massive (M > 2.5 M©) main-sequence (MS) 
stars of spectral type OB harbour a fully developed convec¬ 
tive core below a radiative envelope during their core hydrogen¬ 
burning phase. As their MS evolution progresses, the decrease 
in the central hydrogen content Xc implies that the convective 
core shrinks and it hence leaves behind a chemical gradient 
of CNO-processed material at the interface between the inner 
convective core and the outer radiative envelope. The extent of 
this interface zone and its mixing properties remain largely un¬ 
known. Asteroseismology in principle offers the opportunity to 
probe this zone. 

Stars are multi-component fluids, consisting of chemical 
species i with a mass fraction denoted as Xi. The spatial abun- 
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dance distribution X/(r, t) does not stay constant during the evo¬ 
lution. In burning regions, thermonuclear reactions provide a 
balance between the depletion of the elements i that are burnt 
and the production of heavier elements, while the chemical mix¬ 
ture stays intact in the envelope. In the absence of rotation and 
in a diffusive approximation, the spatial and temporal change 
of the individual mass fractions Xi are governed by the follow¬ 
ing transport equation in the Lagrangian description (e.g. Heger 
et al. 2000): 


dt 


dt), 


+ 


_d_ 

dm 


{Anr^pf Dmix 



( 1 ) 


where the first term accounts for the changes induced by nuclear 
burning, while the second term accounts for the mixing of ele¬ 
ment i according to its spatial gradient dXildm. In this equation, 
Drnix is the effective diffusive mixing coefficient, which has a 
unit of area over time (e.g. Maeder 2009). Advective phenom¬ 
ena such as meridional circulation are excluded in this formu¬ 
lation. Equation (1) is solved subject to boundary conditions for 
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the abundance gradients in the core and at the surface, that is, 
(dXildm)coce = idXildm)sadiice = 0. 

A collection of mechanisms contributes to the local mixing 
in the stellar interior. Convection is the dominant one, because it 
homogenises convective zones on a dynamical timescale, which 
is very short compared to the nuclear timescale. Mixing in the ra¬ 
diative interior occurs on a much longer timescale. A handful of 
mixing mechanisms that can operate in the stellar envelope have 
been proposed, such as core convective overshooting (Saslaw & 
Schwarzschild 1965; Roxburgh 1965; Shaviv & Salpeter 1973), 
semi-convective mixing (Schwarzschild & Harm 1958; Kato 
1966; Langer et al. 1983), rotational mixing (Weiss et al. 1988; 
Charbonnel et al. 1992; Chaboyer 1994), shear-induced mixing 
(Jeans 1928; Zahn 1992), and mixing induced by internal mag¬ 
netic fields (Heger et al. 2005). Following a multivariate analysis 
of the seismic, magnetic, and rotational properties of a sample of 
slowly rotating Galactic OB-type stars, Aerts et al. (2014) found 
the observational evidence for pulsationally induced mixing to 
be stronger than for rotational mixing. 

In the computation of stellar evolution models, most of these 
mixing mechanisms are effectively defined in a diffusion approx¬ 
imation and described by a diffusion coefficient (e.g. Heger et al. 
2000). In the absence of quantitative observational constraints 
on the individual contributions for each of the mixing mecha¬ 
nisms, their net effect ZJrnix is either linearly summed up (Heger 
et al. 2000, 2005) or their interaction is considered (Maeder et al. 
2013; Ding & Li 2014). In this paper, D^ix in radiative interior is 
a linear superposition of individual mixing processes. Any local 
modification to the shape of the mixing profile Dniix(^) propa¬ 
gates into the stellar structure through Eq. (1) and subsequently 
influences the shape of the Brunt-Vaisala frequency, which in¬ 
corporates the gradient of the mean molecular weight. Given the 
sensitivity of high-order g-modes to the detailed shape of the 
Brunt-Vaisala frequency, their frequencies are infiuenced by the 
local shape of Dniix(r). For this reason, and thanks to asteroseis- 
mology based on the high-precision Kepler space-based photom¬ 
etry, we are able to investigate whether observational constraints 
on Dmix are within reach. 

Before high-precision space photometry was available, 
Miglio et al. (2008) have already shown how varying D^ix infiu- 
ences the morphology of the period spacings of gravity modes 
in MS stars with convective cores. Degroote et al. (2010) ap¬ 
plied this approach to model the detected period spacing of the 
slowly rotating pulsating B3V star HD 50230 observed with the 
CoRoT mission for five months. They assumed the detected se¬ 
quence of consecutive modes belongs to the dipole {f = 1) series. 
They concluded that HD 50230 has a mass between 7 and 8 M©, 
has consumed some 60% of its initial hydrogen, and required a 
value of logDmix between 3.48 and 4.30 cm^ s“^ to explain the 
small periodic deviation of 240 s from a constant period spac¬ 
ing of 9418 s for an assumed value of 0.02 for the metallicity. 
However, this result was criticised by Szewczuk et al. (2014), 
who suggested that some of the low-amplitude peaks in the se¬ 
ries might be due to modes of different degree. 

KIC 10526294 was recently discovered to be a very slowly 
rotating and slowly pulsating B (SPB) star from Kepler photom¬ 
etry (Papics et al. 2014, hereafter P14), exhibiting 19 rotationally 
split dipole gravity modes. Papics and collaborators performed 
forward seismic modelling and considered core overshooting 
in the diffusive exponentially decaying prescription defined in 
Freytag et al. (1996) and Herwig (2000). We designate the free 
parameter of this prescription as /ov (and discuss it in Eq. (3) be¬ 
low). But, they only succeeded to put an upper limit on this free 
parameter: /ov < 0.015. 


In contrast to the case of HD 50230, there is no doubt about 
the identification of the degree of the modes for KIC 10526294. 
Hence, we here revisit the forward modelling of this SPB 
(Sect. 2) and introduce two improvements. The first one con¬ 
cerns the inclusion of realistic frequency uncertainties, to com¬ 
pute reduced values as a means of the model selection cri¬ 
terion. P14 only adopted a first rough goodness-of-fit procedure 
for the model selection, denoted here as;^p^ 4 . This goodness-of- 
fit was based on the conservative Rayleigh limit of 0.000685 d“^ 
for the frequency error of all 19 detected modes, with the argu¬ 
ment that the uncertainties on theoretical predictions of oscilla¬ 
tion frequencies of pulsating B stars are typically about 0.001 d“^ 
(Briquet et al. 2007). That was adequate for their modelling, be¬ 
cause they only compared models with one set of input physics 
and with one fixed number of degrees of freedom. Since we 
here compare models with and without extra diffusive mixing, 
and thus with a different number of degrees of freedom, we 
must use a more appropriate statistical description to perform the 
model selection. This requires the determination of individual 
frequency uncertainties, as explained in Sect. 2. The second im¬ 
provement concerns a more detailed model comparison, where 
we consider different elements of the input physics, such as 
opacities, chemical mixtures and two different descriptions for 
the core overshooting. Moreover, we add models based on the 
inclusion of extra diffusive mixing outside the convective and 
overshooting zones, described by a coefficient Dmix- 

We introduce the target star in Sect. 2 and explain how we 
took into account its seismic data in Sect. 3. The physical in¬ 
gredients of the evolutionary models computed with the MESA 
code (Paxton et al. 2011, 2013, version 5548) are explained in 
Sect. 4. The adiabatic zonal dipole frequencies for each input 
model along the evolutionary tracks were computed with the 
GYRE pulsation code (Townsend & Teitler 2013, version 3.0). 
In Sect. 5, we explain the model selection based on a reduced 
method. The results of our grid search to select the best mod¬ 
els to explain the seismic data of KIC 10526294 are reported in 
Sect. 5. Our consideration of three metal mixtures and both OP 
and OPAL opacities is exploited in terms of the mode stability 
for the best models in Sect. 6. Finally, we discuss our results in 
Sect. 7. 

2. The showcase of KIC 10526294 

P14 identified KIC 10526294 as a new SPB star (e.g. Waelkens 
1991; De Cat & Aerts 2002, for a definition) from its Kepler 
light curve, which has a duration of AT = 1460 d. Its effective 
temperature, surface gravity, and metallicity were determined by 
P14 from high-resolution spectroscopy: Teff = 11 500 + 500 K, 
logg = 4.1 + 0.2 dex, and Z = 0.016^^^. P14 performed a full 
study of the Kepler photometry and identified 19 consecutive 
dipole i = I g-mode triplets. As mentioned above, we consider 
these 19 zonal-mode frequencies with their appropriate errors. 
The latter were determined following the methodology discussed 
in detail in Degroote et al. (2009), which is based on the the¬ 
ory of time-series analysis of correlated data by Schwarzenberg- 
Czerny (1991). In practice this implies taking the formal errors 
of the non-linear least-squares fit to the light curve and to correct 
them for the signal-to-noise ratio, sampling, and correlated na¬ 
ture of the data. Applying the procedure discussed in Degroote 
et al. (2009) to KIC 10526294 results in a correction factor of 
3.0 to be applied to the formal errors listed in TableB.l in PI4. 
For convenience, we repeat the 19 zonal dipole frequencies, their 
periods and period spacings and their corresponding Icr errors in 
Table 1. 
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Table 1. Frequencies periods period spacings APf^^\ and 
Icr errors of the central peaks (m = 0) of the dipole = 1) triplets 
observed in KIC 10526294 as determined from Table B.l in P14. 


Id 

' [d-'] 

/i 0.472220+0.000057 
/2 0.486192+0.000057 
/s 0.500926+0.000036 
/4 0.517303+0.000051 
/5 0.533426+0.000018 
fe 0.552608+0.000006 
fi 0.571964+0.000012 
/s 0.593598+0.000039 
fg 0.615472+0.000018 
/lo 0.641202+0.000069 
fn 0.670600+0.000057 
/i2 0.701246+0.000075 
/i3 0.734708+0.000009 
fu 0.772399+0.000045 
/i5 0.812940+0.000009 
fie 0.856351+0.000030 
/i7 0.902834+0.000069 
/i8 0.954107+0.000030 
/i9 1.013415+0.000036 


+ 6 ,- 

I ^ I ^ 

[d] [s] 

2.11766+0.000256 — 

2.05680+0.000241 5258+30 

1.99630+0.000143 5227+24 

1.93310+0.000191 5460+21 

1.87467+0.000063 5048+17 
1.80960+0.000020 5622+6 

1.74836+0.000037 5291+4 

1.68464+0.000111 5505+10 

1.62477+0.000048 5173+10 

1.55957+0.000168 5633+15 

1.49120+0.000127 5907+18 

1.42603+0.000153 5631 + 17 

1.36108+0.000017 5612+13 

1.29467+0.000075 5738+7 

1.23010+0.000014 5578+7 

1.16775+0.000041 5388+4 

1.10762+0.000085 5195+8 

1.04810+0.000033 5143+8 

0.98676+0.000035 5300+4 


As discussed in PI4, the average rotational frequency split¬ 
ting for the m = +1 components of the dipole modes points to a 
rotation period of ~ 188 d when averaged over the entire depth of 
the star, implying that KIC 10526294 is an ultra-slow rotator. Ac¬ 
cording to the survey of galactic B stars by Huang et al. (2010), 
the probability of finding a very slowly rotating late B-type star 
(2 < M/Mo < 4) is low (their Fig. 7a). Thus, KIC 10526294 is 
optimally suited to study mixing mechanisms for a case where 
the effect of rotation is negligible. 

Based on the initial asteroseismic modelling presented in 
P14, the mass of KIC 10526294 is roughly 3.2 Mq, and the star is 
situated very close to the ZAMS. P14 assumed convective core 
overshooting and semi-convective mixing (according to the pre¬ 
scription by Danger et al. 1983, with a^c = 10“^). They found 
their best seismic models (their Table 6) to have a core over¬ 
shooting /ov in the range from zero to 0.015. In what follows, 
we start from the same grid of models as in P14 but consider ad¬ 
ditional extra mixing, various chemical mixtures, different opac¬ 
ity tables and also a step function for the core overshooting. For 
all these cases, we study how the goodness-of-fit behaves and 
perform model comparisons. 


3. Period spacings of the gravity modes 

The probing power of asteroseismology of SPB stars emerges 
from (a) the sensitivity of the g-mode frequencies to the physics 
of the stellar interior in general and to the propagation cavities in 
particular, (b) the relative eigenfunction variation across the star, 
and (c) mode trapping in the overshooting zone (e.g., Dziem- 
bowski & Pamyatnykh 1991; Miglio et al. 2008). At a fixed 
age (associated with the hydrogen mass fraction in the core Xc), 
some modes have a node near the chemically inhomogeneous 
zones and are able to probe the overshooting layer. 


mnn 



Period [day] 


Fig. 1. Top. The period spacing AP for the 19 dipole gravity modes 
observed in the Kepler SPB star KIC 10526294. Bottom. The Icr uncer¬ 
tainty region around each frequency as given in Table 1. 


The asymptotic period spacing for high-order low-degree 
gravity modes for a non-rotating star is given by 


~ Pn+\,i Pn,£ ~ 


+ 1) \ J/?cc ^ / 


( 2 ) 


where is the stellar radius, Rcc indicates the boundary of the 
convective core, and Pn,t and Pn+\,t are periods of two consecu¬ 
tive modes in radial order of the same degree I (Tassoul 1980). In 
practice, the spatial integration is only carried out in the region 
where the g-modes propagate. Since we only consider consecu¬ 
tive dipole modes, we drop the n and i subscripts from APnj in 
the rest of the paper. 

Fig. la presents the sequence of 18 period spacings for the 
19 consecutive dipole g-modes detected in KIC 10526294. They 
exhibit a small but clear deviation from the asymptotic (i.e. con¬ 
stant) period spacing relation in Eq. (2) and highlight the pres¬ 
ence of chemical inhomogeneities that are left behind by the 
shrinking convective core (Miglio et al. 2008). In Fig. lb, the 
frequency uncertainties ct/ of each frequency are presented 
as a function of the mode periods. 

During the MS phase, the core gradually shrinks and its den¬ 
sity and pressure increase. Thus, the integrand in the right-hand 
side of Eq. (2) increases, and the asymptotic value AP as well 
as the structure of the period spacing sequence evolve. The for¬ 
mer depends on the size of the fully mixed core, while the latter 
depends on the detailed treatment of mixing in the radiative en¬ 
velope on top of the shrinking core (Miglio et al. 2008). Fig. 2 
shows the evolution of the asymptotic AP value computed from 
the right-hand side of Eq. (2) for three models with masses of 
2.5 Mo, 3.2 Mq, and 3.8 Mq. The observed period spacing range 
for KIC 10526294 from Table 1 is highlighted in grey. In prac¬ 
tice, we chose the median frequency, /, from Table 1 and carried 
out the integration in the region where (with S£ 

the Lamb frequency). Because the frequency range of /i to fig 
is limited, the choice of / does not alter the result of the inte¬ 
gration, as we have verified. The asymptotic period spacing for 
models without overshooting (/ov = 0.00, solid black lines) are 
mostly below those for models with overshooting (/ov = 0.03, 
blue dashed lines). Varying the initial metallicity Zini introduces 
similar changes, which are not shown here for clarity. 

Because several combinations of the stellar parameters can 
reproduce the same asymptotic AP value, it is essential to fit in¬ 
dividual frequencies in addition to matching the asymptotic AP 
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Core Hydrogen 

Fig. 2. Evolution of the asymptotic period spacing for dipole modes 
computed from the right-hand side of Eq. (2) versus for models with 
masses 2.5 M©, 3.2 M©, and 3.8 M©, each for Zm = 0.020. The solid 
lines correspond to models without overshooting, while the dashed lines 
correspond to models including core overshooting with a value of /ov = 
0.03. The thin grey bar highlights the measured period range of 5048 to 
5907 s for KIC 10526294 from Table 1; see also Eig. 1. All models were 
computed using OPAL opacity tables, and the Nieva & Przybilla (2012) 
and Przybilla et al. (2013) mixture. 


when performing seismic modelling. Based only on Fig. 2, both 
models at the low and high side of the mass range [3,4] M© are 
able to explain the observed AP for KIC 10526294. Additional 
tuning and model selection requires scanning dense asteroseis- 
mic model grids that have a sufficiently broad range in mass M, 
but are also sufficiently fine in terms of all the other input pa¬ 
rameters Aini, Zini, Ac, logDmix, and /ov, to obtain a meaningful 
result, according to the precision of the identified frequencies. 


4. Extensive evolutionary and asteroseismic 
grids 

P14 presented a dense evolutionary and asteroseismic grid of 
models dedicated to KIC 10526294; they considered relatively 
young models (Ac > 0.40) in the range of 3.00 to 3.40 M©. 
Here, we followed the same approach, but we also included extra 
diffusive mixing beyond the convective and overshooting zones. 
Moreover, we varied several input parameters to study if that im¬ 
proves the quality of the frequency fitting compared to the one 
achieved by P14. We present eight evolutionary and asteroseis¬ 
mic grids in which some parameters are varied and some are kept 
fixed to keep the CPU requirements manageable. 

The physical ingredients of the grids are the following. We 
used the open-source MESA^ code (Paxton et al. 2011, 2013) to 
calculate evolutionary tracks and store equilibrium models along 
each track for every 0.001 drop in Ac. This tiny grid step in the 
central hydrogen fraction was taken to ensure a small enough 
change in the frequency values of the modes as we follow the 
stellar evolution along an evolutionary track, relying on the nu¬ 
merical accuracy achieved by the MESA code coupled to the 
pulsation code we used (as discussed below), cf. Fig. 13 in P14 
and its discussion. 

We adopted the Ledoux convection criterion with a high 
semi-convective mixing coefficient given by a^c = 10“^ in 
the prescription by Langer et al. (1983). For this value, the 

^ MESA can be downloaded from the project Web site 
http://mesa.sourceforge.net 


Ledoux criterion is equivalent to the Schwarzschild criterion 
and ensures that Vj-ad = Vad on the convective side of the core, 
following Gabriel et al. (2014). The mixing length parameter 
was set to ctmlt = 1-8, while we fixed the initial composi¬ 
tion to the Galactic standard given by Nieva & Przybilla (2012) 
and Przybilla et al. (2013, hereafter NP12): (Aini, Uni, Zini) = 
(0.710,0.276,0.014). We used the OPAL Type 1 (Rogers & Nay- 
fonov 2002) opacity tables adapted to this mixture; these tables 
and the MESA inlists are available for download as online mate¬ 
rial (see Appendix C). 

We adopted an exponentially decaying prescription for over¬ 
shooting on top of the convective core in a diffusive approxima¬ 
tion, following Freytag et al. (1996) and Herwig (2000): 


^ov(^) — ^conv 



(3) 


where L>conv is the convective mixing coefficient, z = r - R^c is 
the distance above the boundary of the convective core Rcc into 
the radiative region, and Hp is the local pressure scale height. 
The parameter /ov governs the width of the overshooting layer, 
which in the framework of the mixing length theory of Bohm- 
Vitense (1958), is not constrained from first principles. We there¬ 
fore treated it here as a free parameter and varied it in a broad 
range from 0.0 to 0.03 in all our grids. 

An alternative for core overshooting is to use a step-function 
prescription, which is commonly used in the literature. In this 
prescription, the overshooting extends over a distance Jov = 
aoyHp, beyond the convective core boundary. The material is 
treated as fully mixed in the overshoot range: 


Doy 


= D 


conv? 


Rqc — — Rcc ^ov 


(4) 


where Rcc is the radial position of the boundary of the convective 
core, Jov is the width of the overshooting zone and L>conv is the 
diffusion coefficient in the fully mixed convective core. Follow¬ 
ing the MESA implementation of core overshooting according 
to a step function, we tested explicitly that our seismic models 
are insensitive to the depth in the convective core at which the 
value of Dconv is picked up. 

Our eight extensive evolutionary and asteroseismic grids are 
composed and named as follows: 

- Basic Grid: We investigated whether more massive and 
more evolved models may reproduce the observed, nearly 
fiat, period spacing pattern of KIC 10526294 as well. The 
Basic Grid is essentially the extension of the grid in P14 to¬ 
wards higher masses and lower Ac. 

- Composition Grid: We calculated a grid varying Aini and 
Zini for the NP12 mixture to investigate if this improves the 
frequency fit. 

- Mixing Grid: Various extra mixing mechanisms can operate 
and lead to smoothing of the composition profile outside the 
core. Unfortunately, we have no solid quantitative measures 
of the amount of extra mixing for B stars. Hence we com¬ 
puted a grid with models including extra diffusive mixing 
throughout the star for a wide range of log Dmix values. 

- Evolved Grid: A large amount of extra mixing combined 
with a low Ac may provide a situation for which the AP 
profile matches the observations for higher mass models. In 
Fig. 2, higher mass models cross the grey band at lower Ac. 
To ensure our age constraint is reliable, we calculated a grid 
for M > 3.30 M© incorporating extra mixing; we only con¬ 
sidered the models whose mean AP falls inside the observed 
range in Fig. 2. Thus, the starting and ending Ac values for 
each track are fiexible. This is the meaning of “—” in Ta¬ 
ble 2. 
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Fig. 3. Comparison between observed periods (vertical lines) and those of six models with different mixtures and opacity tables. All parameters 
of the models are identical, except their mixture and opacities. The widths of the vertical lines are larger than the measured pulsation period 
uncertainties. 


- Fine Grid: To resolve the parameter space of mass, over¬ 
shooting, and extra mixing, we increased the resolution for 
these parameters around the best model found from the Mix¬ 
ing Grid. Thus, we varied M e [3.18, 3.27] M© in steps 
of 0.01 Mq, /ov € [0.016,0.019] in steps of 0.001, and 
logDmix ^ [0.25, 2.75] in steps of 0.25 dex. 

- Metallicity Grid: Although the average metallicity of 
Galactic B stars in the solar neighbourhood is found to be 
Z = 0.014 + 0.002 (Nieva & Przybilla 2012), individual B 
stars may still show a strong deviation from this average. 
Enhanced metallicity increases the height of the iron opacity 
bump and ensures better mode excitation. Moreover, as we 
demonstrate below in Fig. 5, the initial metallicity and mass 
exhibit a strong anti-correlation. To explore the possibility 
of our target having much higher metal abundance than the 
cosmic standard, we computed a grid identical to the Mixing 
Grid and Fine Grid, but with lower initial mass and higher 
initial metallicity up to Zini = 0.030. 

- Mixture Grids: The theoretical pulsation frequencies de¬ 
pend on the initial metal mixture and the adopted opacity 
tables. To illustrate this. Fig. 3 compares the observed peri¬ 
ods (vertical grey lines) with those of six models with iden¬ 
tical model parameters, except that they have different mix¬ 
tures and opacities. The adopted mixtures for the grid are 
NP12, the solar mixture reported by Asplund et al. (2009) 
(hereafter A09) and the Asplund et al. (2005) mixture with 
a Ne enhancement based on Cunha et al. (2006) (hereafter 
A05-fNe). Fig. 3 shows a clear difference in the frequency 
predictions, particularly for the longer-period modes. The 
relative difference in periods ranges from 0.07% to 0.14% 
from left to right in the figure. Changing the mixture in the 
models infiuences fitting high-precision frequencies of well- 
identified pulsation modes. In this respect, the asteroseismic 
data of KIC 10526294 bring us to the level where the depen¬ 
dence of the models on the assumed mixture can be tested 
for SPB stars, thanks to the identification of its dipole mode 
triplets. 

We calculated grids using both OP (Seaton 1987) and OPAF 
(Rogers & Nayfonov 2002) opacity tables for each of the 
NP12, A09, and A05-i-Ne mixtures. The parameter range for 
each of these sub-grids is identical. Once we found the best 
model in each sub-grid, we repeated the computations within 
a smaller range around the best model until the step size 
in mass, overshooting, metallicity, and logDmix was 0.01, 
0.001, 0.001 and 0.25 dex, respectively. 

- Step Function Overshoot Grid: To examine which of the 
two prescriptions for the core overshooting, exponentially 


decaying as in Eq. (3) versus step function as in Eq. (4), pro¬ 
vides a better fit to the observed data, a stand-alone grid of 
models for a step function overshoot was computed around 
the best model from the Fine Grid, with the only change of 
using the step function instead of the exponential prescrip¬ 
tion for a range of masses and metallicities. 

Table 2 gives the full parameter range in each of these grids. 
In the Mixing Grid and High Mass Grid, we varied the minimum 
diffusive mixing coefficient Z9mix- For the Basic Grid and Com¬ 
position Grid, we set Dmix = 0. In all grids except the Step Over¬ 
shoot Grid, we used the exponential prescription as in Eq. (3). In 
total, our grids comprise 27 058 evolutionary tracks and around 
4.3 million models that serve as input for the GYRE pulsation 
computations. 

For each input model, we used the state-of-the-art GYRE 
code (Townsend & Teitler 2013, version 3.0) to solve for the 
linear adiabatic oscillation frequencies of dipole ^ = 1, m = 
0 modes in the frequency range of the detected g-modes of 
KIC 10526294. The input inlist for the GYRE computations is 
available for download as online material; see Appendix C. 


5. Model selection 

5.1. Methodology 

We considered a reduced minimisation scheme to select 
the best-fitting models after comparing the observed frequencies 
with their model counterparts : 


2 

Yred 


1 /f(obs) Ath)x2 

i y G/ ~ h I 

Nf - n ^ cfi 

j i=\ I 


(5) 


where Nf = 19 is the number of consecutive g-modes and ct/ 
are the theoretical frequency uncertainties listed in Table 1. The 
number of degrees of freedom, n in Eq. (5), represents the num¬ 
ber of independent model parameters. In all the cases, n is either 
4 or 5, as listed in Table 2. This procedure is a simplified version 
of the goodness-of-fit adopted by PI4, who also included the pe¬ 
riod spacing values in the;^^. Here, we prefered to forego this 
because demanding an optimal fit to the observed frequencies 
individually will automatically also deliver an optimal fit to the 
observed period spacing values. Moreover, we wished to com¬ 
pare the fit quality of models with a different number of degrees 
of freedom (^ = 4 or 5, cf. Table 2) for the best model selection 
and this is better justified statistically when using independent 
measurements in the computation of x^^^- 
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Table 2. Eight model grids and their parameter range dedicated to the 
asteroseismic modelling of KIC 10526294. N is the number of values 
for each parameter. cTsc is fixed to 10"^ in all models. Dmix is in cm^ s"^ w 
is the number of degrees of freedom in Eq. (5). Eor explanations of in 
the Evolved Grid, we refer to the text. The Mixture Grid is repeated for 
OP and OPAL opacity tables and for NP12, A09 and A05+Ne mixtures. 


Parameter 

Erom 

To 

Step 

N 

Basic Grid {n=4): 

Ini. mass: Mm [M©] 

3.00 

3.95 

0.05 

25 

Overshooting: /ov 

0.000 

0.030 

0.003 

11 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.010 

0.020 

0.001 

11 

Centre hydrogen: Xc 

0.700 

0.0 

0.001 

>701 

Composition Grid {n=5)\ 

Ini. mass: Mm [M©] 

3.15 

3.25 

0.05 

3 

Overshooting: /ov 

0.000 

0.030 

0.003 

11 

Ini. hydrogen: Xm 

0.68 

0.72 

0.01 

5 

Ini. metallicity: Zm 

0.010 

0.020 

0.001 

11 

Centre hydrogen: Xc 

0.700 

0.600 

0.001 

101 

Mixing Grid {n=5)\ 

Ini. mass: Mm [M©] 

3.10 

3.30 

0.05 

5 

Overshooting: /ov 

0.000 

0.030 

0.003 

11 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.010 

0.020 

0.001 

11 

Centre hydrogen: Xc 

0.700 

0.600 

0.001 

101 

Extra Mixing: log Dmx 

2 

6 

1 

5 

Evolved Grid {n=5)\ 

Ini. mass: Mm [M©] 

3.30 

4.00 

0.05 

15 

Overshooting: /ov 

0.000 

0.030 

0.006 

6 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.010 

0.020 

0.002 

6 

Centre hydrogen: Xc 

— 

— 

0.001 

— 

Extra Mixing: log Dmx 

2 

6 

1 

5 

Fine Grid {n=5)\ 

Ini. mass: Mm [M©] 

3.18 

3.27 

0.01 

10 

Overshooting: /ov 

0.000 

0.030 

0.003 

11 


0.016 

0.019 

0.001 

3 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.010 

0.020 

0.001 

11 

Centre hydrogen: Xc 

0.700 

0.600 

0.001 

101 

Extra Mixing: log Dmx 

0.25 

2.75 

0.25 

10 

metallicity Grid {n=5)\ 

Ini. mass: Mm [M©] 

3.00 

3.20 

0.05 

5 


3.12 

3.18 

0.01 

7 

Overshooting: /ov 

0.012 

0.024 

0.003 

5 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.021 

0.030 

0.001 

11 

Centre hydrogen: Xc 

0.700 

0.600 

0.001 

101 

Extra Mixing: log Dmx 

1.00 

2.00 

0.25 

5 

Mixture Grids {n=5)\ 

Ini. mass: Mm [M©] 

3.17 

3.29 

0.02 

7 

Overshooting: /ov 

0.012 

0.024 

0.003 

5 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.010 

0.018 

0.002 

5 

Centre hydrogen: Xc 

0.700 

0.600 

0.001 

101 

Extra Mixing: log Dmx 

1.50 

2.50 

0.25 

5 

Step Overshoot Grid {n=5)\ 

Ini. mass: Mm [M©] 

3.18 

3.27 

0.01 

10 

Overshooting: Ofov 

0.16 

0.23 

0.01 

8 

Ini. hydrogen: Xm 

0.71 

0.71 

0.01 

1 

Ini. metallicity: Zm 

0.010 

0.020 

0.001 

10 

Centre hydrogen: Xc 

0.700 

0.600 

0.001 

101 

Extra Mixing: log Dmx 

1.50 

2.25 

0.25 

4 


Thex^Qd in Eq. (5) gauge how well the model frequen¬ 

cies match the observed ones, where a value of 1 would mean 


that the input physics of the models is correct up to the level of 
the frequency precision. We are very far from this situation in the 
case of stars of spectral type O, B, or A. Indeed, for such mas¬ 
sive stars, we are in a different situation compared to solar-like 
stars with stochastically excited oscillations, for which scaling 
relations can be applied to the Kepler data (Chaplin et al. 2014) 
and where extrapolations from the solar model are meaningful. 
In these straightforward cases, one achieves seismic modelling 
with typical values of x^^d ~ Kepler data of 

G- or F-type stars on the main sequence, is achieved with a grid- 
based approach like the one we adopt here (Metcalfe et al. 2014). 
For the best modelled subgiants with solar-like oscillations, 
values between 100 to 2000 were obtained (e.g. Deheuvels et al. 
2014). 

Seismic modelling of massive stars with heat-driven modes 
has not been done so far using the level of precision of the mea¬ 
sured frequencies, cf. the discussion in Sect. 1 concerning the 
use of the Rayleigh limit as rough estimate of the measured fre¬ 
quency precision. Theoretical frequency predictions from mod¬ 
els have typically only been computed up to 0.001 d“\ and the 
models were selected by visual inspection of the frequency spec¬ 
tra (e.g. Pamyatnykh et al. 2004) or by a;^^-type goodness-of-fit 
function ignoring the measured frequency precisions (e.g. Bri¬ 
quet et al. 2007). This is appropriate as a procedure when only a 
handful of identified mode frequencies is available. In this way, 
stellar masses, radii, core overshooting values, and ages have 
been deduced (e.g. Aerts 2015, for an overview). 

5.2. Candidate models 

Tables lists the parameters of the best-fit models for each of 
the grids described in Table 2. The first and second columns as¬ 
sign a number to each model and specify their grid of origin, 
respectively. The next eight columns specify the adopted input 
parameters for the best models. The last two columns give the 
X^^^ (Eq. 5) and Xpu values, respectively. The latter column is 
given to provide a sensible comparison with the results obtained 
from our previous modelling as discussed in P14 (their Table 6, 
Col. 11), where the input physics was kept fixed. 

Model 1 is essentially the replication of the best model found 
in P14 from the Basic Grid, as a verification of the consis¬ 
tency between our previous and current study. Model 2 is taken 
from the Composition Grid. Comparison between Model 1 and 
Model 2 shows that, for the specific case of KIC 10526294, 
changing the initial composition Aini and Fini deteriorates the 
frequency fittings. Thus, the Galactic standard composition of 
B stars by Nieva & Przybilla (2012) is the most appropri¬ 
ate one to take for this target. Models is the best of its kind 
from the Evolved Grid, but its x^^d significantly higher than 
that of Model 1. Therefore, we exclude the possibility that 
KIC 10526294 is an evolved star. 

Model 4 is taken from the Mixing Grid and the Fine Grid. 
It has the lowest x^^d is t)^st seismic model for 

KIC 10526294. It outperforms all models without diffusive mix¬ 
ing and is roughly 11 times better in fit quality than Model 1. 
Thus, we conclude that the frequency fitting for KIC 10526294 
requires extra diffusive mixing with a value of logD^ix = 
1.75 cm^s“^ This is a similar conclusion as was reached for 
HD 50230 (Degroote et al. 2010), only this time we were able 
to quantify D^ix much more precisely and we considered vari¬ 
ous options for the metal mixture and opacities, while Z = 0.02 
was fixed for HD 50230 and the degrees of its modes were not 
identified observationally. We find a value for the diffusive mix¬ 
ing coefficient of KIC 10526294 that is two orders of magnitude 
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lower than the one for HD 50230. The two major differences be¬ 
tween these stars are their mass and their Xc, both being ultra- 
slow rotators. It would be worth revisiting the seismic modelling 
of HD 50230 now that we are aware that it is a member of a spec¬ 
troscopic binary and that p-modes have been found in addition to 
g-modes (Degroote et al. 2012); these two facts were not taken 
into account in the modelling efforts made by Degroote et al. 
( 2010 ). 

Models 5 to 9 originate from the Mixture sub-grids. They 
provide a test of the influence of different opacities and mixtures 
on the best values for /ov and logD^ix- Regardless of their 
Tred all these models converge to very similar values for 

Xc=0.63 to 0.64, /ov=0.017 to 0.018, and logDniix=1.75 to 2.00. 
However, the mass varies somewhat and the metallicity varies 
widely over the scanned range. The fact that X^ /ov and log D^ix 
hardly depend on the mass, metallicity, mixture, and opacities 
manifests that a solid seismic constraint is placed on these values 
for the speciflc case of KIC 10526294. The discretisation of the 
grid parameters is also sufficiently optimal that the best values 
found are close to one another. 

Models 5 and 8 have the A09 mixture and provide the fourth 
and third best fits to the frequencies, respectively. Although they 
are built from different opacities, the fit quality of these two 
models is indistinguishable in a statistical sense. Both are of 
lower fit quality than Model 4, given that the difference between 
their of Model 4 is far above 3.84 (which would 

correspond to a p-value of 0.05 in the case of nested models). 
Furthermore, Model 4 and Model 7 have identical composition 
and mixture, but were calculated using OPAL and OP opacity ta¬ 
bles, respectively. The parameters of these two models are very 
close in fit quality, although the latter has the worse Mod¬ 
els 6 and 9, with the A05-fNe mixture, also provide reasonably 
good fits to the frequencies. Therefore, changing the mixture, at 
least in this case, influences the adiabatic frequencies more than 
swapping between OP and OPAL opacity tables. 

The x^Qd Model 10 is more than twice as high as the one 
of the best model. Hence, we come to the important conclusion 
that, for this slowly rotating SPB star the exponentially decay¬ 
ing prescription for core overshooting gives a better fit to the 
seismic data than a step function. A more reliable distinction 
between the two prescriptions for massive pulsators in general 
requires a reasonably large ensemble of B-type pulsators to be 
modelled with the two prescriptions. A re-analysis of the sam¬ 
ple in Aerts (2015) could be highly instructive in this regard. 
Previously, Dziembowski & Pamyatnykh (2008) already consid¬ 
ered an overshooting prescription based on two free parameters 
rather than a simple step function formulation in their modelling 
of the jSCep stars vEri and 12 Lac, but these two pulsators did 
not allow Arm conclusions on the shape of the overshooting. 

Model 11 is selected from the Metallicity Grid, and is ranked 
as the second best asteroseismic model for our target. Even 
though it has a different initial mass and metallicity, the over¬ 
shooting and extra mixing parameters perfectly agree with those 
of Model 4 to Model 9. Because Model 11 has 100% higher 
metal abundance than Model 4, it has a distinct stability prop¬ 
erty that is addressed in Sect. 6. Based on our;^^^^ minimisation. 
Model 4, and to a lesser extent Models 11 and 8, best represent 
the asteroseismic data of KIC 10526294. In the following sec¬ 
tions, we use these three models for illustrative purposes. 

5 . 3 . Comparing best models and observations 

Fig. 4a compares the observed period spacing (grey symbols) 
versus mode periods with those of Model 4 (fllled circles) and 


Model 1 of P14 (empty squares), respectively. The resulting pe¬ 
riod spacing pattern of Model 4 clearly follows the observed pat¬ 
tern better than the pattern from Model 1 does. To appreciate how 
well the frequencies from Model 4 approach the observed ones, 
we show in Fig. 4b the frequency deviation Sft between the two 
through 

/=1, •••,19. (6) 

In Fig. 4b, the ordinate range is 40 times larger than in Fig. lb, 
and the Icr frequency uncertainty band appears as a line here. 
The frequencies of Model 1 (empty squares) show stronger de¬ 
viations from zero and the frequencies of Model 4 provide far 
better fits to the observed period series, especially in the lower 
period regime. To appreciate the quality of the model fit, we 
computed the average 6ft for Model 4 and obtained 0.00051 d“/ 
While this is still a factor ~ 30 larger than the observational un¬ 
certainty on the frequencies, it is by far the best seismic model 
constructed to represent the pulsations of a late B-type star up to 
now. 

Additional comparisons of the period spacing fits between 
Models 4 and 11 and between Models 8 and 11, are presented in 
Fig. A.l and Fig. A.2 in the Appendix. We list the global param¬ 
eters of the three best seismic models of KIC 10526294 in Ta¬ 
ble 4. Their structure variables (in GYRE-compatible format) are 
available online; see Appendix C. Even for optically bright B- 
type pulsators, the seismically derived effective temperature and 
surface gravity do not necessarily agree with their spectroscopic 
counterparts because the latter suffer from systematic uncertain¬ 
ties (e.g. Briquet et al. 2007, 2011). In the case of KIC 10526294, 
which is a faint star, the seismic values of the best models all 
agree with the spectroscopic values within 3(T, which we regard 
as good compatibility between seismology and spectroscopy. 

Fig. 5 shows the correlation diagrams between the flve pa¬ 
rameters of the Mixing Grid and the Fine Grid centred around 
the parameters of Model 4, that is, all panels present a subsur¬ 
face in the x^^d where the parameters that are not shown 
on the axes are flxed to the values found for Model 4. This is 
similar to Fig. 11 in PI4, where semi-linear trends between all 
parameters occur. Here, we achieved more stringent constraints 
on some of the parameters, especially X^ /ov, and logD^ix (see 
also Fig.B.l). On the other hand, M and Z are less well con¬ 
strained, although we achieved a relative precision below 4% for 
the mass. A strong linear correlation occurs between them, cf. 
panel (e), which is similar to Fig. 2 in Ausseloos et al. (2004). 
For a full overview of the x^^d distribution for the entire Mixing 
Grid, Fine Grid and Metallicity Grid, we refer to Fig. B.l in the 
Appendix. The existence of correlations between the parameters 
hinders the derivation of appropriate uncertainties on the derived 
parameters of the best model(s), other than the bare minimum 
uncertainty given by the parameter stepsizes. A pragmatic esti¬ 
mate of the overall systematic and statistical uncertainty is pro¬ 
vided by the ranges of the parameters of the best three models 
listed in Table 4. 

Based on Models 4, 8 and 11, we deduce that KIC 10526294 
is a young star (Xc ^ 0.63 to 0.64) of late-B spectral type. 
From the parameters of Models 4 to 9 and Model 11 in Table 3, 
the core overshoot parameter is tightly constrained in the range 
/ov = 0.017 to 0.018. The overshooting zone has enough spatial 
and mass extension to allow partial trapping of some of the pul¬ 
sation modes. This offers a unique opportunity for an in-depth 
study of the properties of the overshooting layer. Additionally, an 
extra diffusive mixing in the radiative part of the star of an am¬ 
plitude logDmix = 1.75 to 2.00 cm^ s“^ is needed to explain the 
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Table 3. Compilation of best models from different grids (Table 2) and their input parameters. These models have the lowest scores within 
each grid. The overshooting parameter is /ov for all models in Col. 9, except for Model 10 where we list cTov Models 4, 8 and 11 represent the best 
three seismic models of KIC 10526294; for each of them, the radial orders of the 19 detected modes range from 14 to 32. The overall best model 
in a statistical sense is indicated in bold. 


# 

Grid 

Opacity 

Mixture 

Mass 

[Mo] 

^ini 

^ini 


> 

o 

log DjTcdx 
[cm^ s“^ 

] 

Api4 

1 

Basic 

OPAL 

NP12 

3.20 

0.71 

0.020 

0.693 

0.000 

— 

18192 

10.9 

2 

Composition 

OPAL 

NP12 

3.15 

0.69 

0.010 

0.662 

0.030 

— 

11402 

23.3 

3 

Evolved 

OPAL 

NP12 

3.50 

0.71 

0.010 

0.254 

0.030 

2 

23 199 

21.8 

4 

Mixing-hFine 

OPAL 

NP12 

3.25 

0.71 

0.014 

0.627 

0.017 

1.75 

1711 

1.42 

5 

Mixture 

OPAL 

A09 

3.18 

0.71 

0.021 

0.640 

0.018 

2.00 

2 500 

3.60 

6 

Mixture 

OPAL 

A05-fNe 

3.20 

0.71 

0.013 

0.641 

0.018 

2.00 

5 905 

4.53 

7 

Mixture 

OP 

NP12 

3.24 

0.71 

0.014 

0.636 

0.017 

1.75 

3 034 

5.23 

8 

Mixture 

OP 

A09 

3.22 

0.71 

0.015 

0.638 

0.018 

2.00 

2 348 

3.23 

9 

Mixture 

OP 

A05-fNe 

3.17 

0.71 

0.012 

0.632 

0.018 

1.75 

3 249 

2.90 

10 

Step Overshoot 

OPAL 

NP12 

3.19 

0.71 

0.019 

0.628 

0.21 

1.75 

3 792 

5.54 

11 

Metallicity 

OPAL 

NP12 

3.13 

0.71 

0.028 

0.628 

0.018 

1.75 

2187 

2.32 


6000 


• Model #4 t 



□ Model #1 1 

5800 


$ KIC 10526294 


u 

q; 

Gh 


5600 - 


5400 


5200 - 


5000 




Fig. 4. (a) Period and period spacings for Model 4 (filled circles) and Model 1 (empty squares) from Table 3. The observed pattern is shown in grey, 
(b) The frequency difference between model frequencies and those detected in the observations Sfi from Eq. (6). Compared to Fig. lb, the ordinate 
is enlarged 40 times; the grey band around zero is the Icr frequency uncertainty range ct/. A similar plot for Model4, Models and Model 11 is 
shown in Figs. A.l and A.2 in the Appendix. 


observed frequencies. The parameters logDmix, /ov, and Xc only 
show a minor dependence on the adopted mixture and choice of 
the opacity tables. 

One may question the success of our forward modelling for 
KIC 10526294, specifically given the high;^^^^ scores of the best 
model(s) in Table 3. We argue that despite such high values, the 


relative frequency deviations between the observed frequencies 
and those of the best seismic model remain below 0.4%. The left 
panels in Fig. 6 show versus mode period for 

Model 4 (top) Model 8 (middle) and Model 11 (bottom), respec¬ 
tively. The strongest deviations occur for modes with periods ex¬ 
ceeding ~ 1.68 days, that is, /i to /g in Table 1. For Model 4, the 
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1.75 2.00 2.25 
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6.0 

5.6 
5.2 
4.8 
4.4 
4.0 

3.6 


3.18 3.21 3.24 3.27 3.18 3.21 3.24 3.27 1.50 1.75 2.00 2.25 

Mass [Mq] Mass [M^] 

Fig. 5. Correlation diagram of the parameters of the Fine Grid around Model 4. Results are based on the Mixing Grid plus the Fine Grid. The 
colour coding is based on log;^^^^^ . To see the improvement we achieved in confining the grid parameter space, compare with Fig. 11 in PI4. In 
each panel, the model parameters that are not on the axis are fixed to those of Model 4. Here, Xc, /ov and log are well constrained. 


Table 4. Global parameters of Model4, Models, and Model 11. The 
free grid parameters and the derived model parameters are separated by 
the horizontal line. Jqv and Mqv are the width and the mass contained in 
the overshooting zone, respectively. 


Parameter 

Model 4 

Model 8 

Modem 

Mass [Mq] 

3.25 

3.22 

3.13 

/ov 

0.017 

0.018 

0.018 


0.627 

0.638 

0.628 

Xini 

0.710 

0.710 

0.710 


0.276 

0.275 

0.262 

Zini 

0.014 

0.015 

0.028 

logDmix [cm^S"'] 

1.75 

2.00 

1.75 

Radius [R©] 

2.215 

2.195 

2.382 

Luminosity [L©] 

128 

111 

78 

Teff [K] 

13 000 

12 650 

11100 

log g (cgs) 

4.259 

4.263 

4.204 

Age [Myr] 

63.0 

61.9 

91.7 

Mec [M©] 

0.672 

0.651 

0.591 

^cc [Rq] 

0.329 

0.326 

0.318 

Mov [M©] 

0.190 

0.195 

0.178 

^OV [Rq] 

0.037 

0.038 

0.037 


first eleven mode frequencies are very similar to the observed 
ones; the deviations gradually grow by increasing mode order. 
Current seismic models of solar-like stars observed with Kepler 
give rise to relative deviations of about 10“"^ between observed 
and modelled frequencies near the frequency of maximum power 
(Metcalfe et al. 2014), but this is only achieved after applying an 
arbitrary surface correction (Kjeldsen et al. 2008). Our results 
are only an order of magnitude worse, but we are treating a star 
very different from the Sun; this is a major achievement. For 
completeness, we point out that we did not find any correlation 
between the mode amplitudes and the relative frequency devia¬ 
tions shown in Fig. 6. 

The right panels in Fig. 6 show histograms of the relative 
frequency deviations between the measurements and theoretical 
predictions around zero. Despite the low-number statistics, the 
deviations for Model 4 are normally distributed, while those of 
the other two models do not show a clear Gaussian shape. For 
the solar-type pulsators, a frequency-dependent correction is ap¬ 
plied to the theoretical frequencies to account for near-surface 
effects (e.g. Kjeldsen et al. 2008). For B-type pulsators, such a 
correction has so far been assumed to be unnecessary because 
these stars have a radiative envelope. If such surface a correction 
were needed for B-type stars, then we would have observed a 
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Period [day] 


Fig. 6. Left panel: Relative deviation between observed and mod¬ 
elled frequencies - y‘(ot's)yy‘(obs) dipole g-modes in 

KIC 10526294. The deviation for all models is below 0.4%. 

Right panel: Histogram distribution of the relative deviations 
around zero. 


monotonically deviating behaviour between the model frequen¬ 
cies and the observed ones, which is clearly not the case. There¬ 
fore, we argue that even if a surface correction term would be 
needed for KIC 10526294, it would be frequency-independent. 

Fig. 7 compares the profile of the Brunt-Vaisala frequency 
N(r), as used in Eq. (2) for four selected models. The inset is 
a zoom-in around the sharp increase in the profile, induced by 
the gradient of the mean molecular weight outside the reced¬ 
ing core. The profiles for Model 4 (black solid). Model 8 (red 
dashed) and Model 11 (blue dotted) are almost identical, be¬ 
cause of their very similar overshooting parameter /ov = 0.017 
to 0.018, except for a slight shift in the position of the peak that 
is due to their different initial masses. The profile for Model 10 
(grey solid) shows a steeper rise around 0.85 M© because of the 
difference between the exponentially decaying and step function 
prescriptions for overshooting. The seismic frequency fitting in¬ 
dicates that the smoother shape of N(r) connected with the expo¬ 
nentially decaying prescription is to be preferred over the steep 
step-function for KIC 10526294. 


6. Mode excitation 

In their Fig. 14, P14 already showed a mismatch between the pre¬ 
dicted and observed excited modes using their best model (i.e. 
Model 1 here). Only eleven modes out of nineteen in Model 1 
are excited. Aside from the forward modelling, we examined if 
our improved modelling remedies that shortcoming. For this, we 
used the non-adiabatic framework of GYRE (version 3.2.2). 

It is well known that current predictions of mode excitation 
through the heat mechanism are not yet sufficiently appropriate 
to explain all the detected and identified modes in B-type stars 
(Dziembowski & Pamyatnykh 2008). For this reason, we exam¬ 
ined the mode stability properties of our best models a posteriori. 

Miglio et al. (2007) already suggested that adopting OP 
opacity tables instead of OPAL might remedy the mode exci¬ 
tation problem. As a result of the significant role of the opac¬ 
ity profile in mode excitation, we first considered its profile in 
three of the best selected models. Fig. 8 compares the Rosse- 
land mean opacity profile k versus temperature log T in Model 4 
(black solid line, using OPAL), Model 8 (red dashed line, using 
OP) and Model 11 (blue dashed-dotted line, using OPAL). The 
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Fig. 7. Bmnt-Vaisala frequency for Model 4 (black solid), Model 8 (red 
dashed), Model 10 (grey solid) and Model 11 (blue dotted) from Ta¬ 
ble 3 and Table 4. The inset is a zoom-in around the peak in the profile 
induced by the gradient of the mean molecular weight outside the con¬ 
vective core. These four models have different masses of the convective 
core Mcc- 


inset is a zoom in of the Fe opacity bump. Obviously, Model 11, 
with its twice as high initial metallicity Zini is more opaque than 
the other two. Moreover, Model 8 is more opaque than Model 4 
despite their nearly identical Zini, because the former model uses 
the OP opacities, in line with Miglio et al. (2007). 

In Fig. 9, we compare the growth rate (i.e. the imaginary part 
of the complex eigenfrequency tuim) as a function of mode period 
for Model 4 (black circles). Model 8 (red squares) and Model 11 
(blue stars), respectively. The vertical lines are the observed peri¬ 
ods. The excited (damped) modes are shown with filled (empty) 
symbols. Although Model 4 reproduces the real part of eigenfre- 
quencies better than Model 11, the latter predicts more excited 
modes than observed. This is not a surprise given the high Zini 
of Model 11, which implies a ~ 70% higher iron opacity (see 
Fig. 8). We point out that an ad-hoc local increase in the height 
of the iron opacity bump by up to ~ 50% (cf. Pamyatnykh et al. 
2004) would allow Model 4 (as well as other models from the 
Mixing Grid and Fine Grid) to have equally good mode excita¬ 
tion properties as Model 11. Model 8 has two additional excited 
modes compared to those of Model 4, despite both having sim¬ 
ilar Zini. This is as expected because OP opacities predict more 
excited modes (Miglio et al. 2007). Additionally, all three mod¬ 
els presented here predict two to four short-period modes to be 
excited but these are not among the detected modes listed in PI4. 

The distinction between Model 4, Model 8 and Model 11 is 
not clear-cut; the former offers the best adiabatic frequency 
match with the observations, while the latter explains the mode 
excitation better at the cost of a somewhat poorer frequency fit¬ 
ting. Additionally, the metallicity value Z^ of Model 4 and of 
Model 8 perfectly agrees with the Galactic standard composition 
of Nieva & Przybilla (2012) for nearby B stars, while Model 11 
is twice as metal-rich as the standard. Given the current lack of a 
higher signal-to-noise spectrum in addition to the one presented 
in P14 upon which we relied, we cannot evaluate this metallic¬ 
ity difference between the best models from spectroscopic abun¬ 
dance analysis. Fortunately, this degeneracy does not affect the 
seismic estimates of the core overshooting and diffusive mixing 
deduced for the star. 
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Fig. 8. Rosseland mean opacity profile in Model 4 (Zini = 0.014, 
OPAL, NP12, black solid line). Model 8 (Zini = 0.015, OP, A09, red 
dashed-dotted line) and Model 11 (Zini = 0.028, OPAL, NP12, blue 
dashed line) versus temperature. The inset is a zoom around the iron- 
bump of opacity. The iron opacity peak for the OP model occurs at 
slightly hotter interior. 



Period [d] 


Fig. 9. Non-adiabatic instability study of Model 4 (black circles). 
Model 8 (red squares), and Model 11 (blue stars) from Table 3. The 
imaginary part of the eigenfrequency dUim is plotted against the 
mode period (vertical lines). Filled (empty) symbols designate excited 
(damped) modes. 


7. Discussion and conciusions 

Forward seismic modelling based on the 19 detected and iden¬ 
tified dipole gravity modes of the B8.3V star KIC 10526294 led 
to the derivation of its mass, radius, and age within the ranges 
M € [3.13,3.25] Mo, R e [2.19,2.38] R©, and age [62,92] Myr, 
respectively. These ranges cannot be refined to tighter constraints 
as long as we do not have a constrained metallicity from spec¬ 
troscopy. In practice, the range of Z e [0.014,0.028] leads to 
good seismic models in this mass range, and the (M,Z) cor¬ 
relation implies a quite large uncertainty on the seismic age. 
However, the seismic data do require the inclusion of a small 
but significant amount of extra diffusive mixing in the radia¬ 
tive envelope of the star, with a value of log Dmix in the range 
[1.75,2.00] cm^ s“\ in addition to an exponentially decaying 
core overshooting with parameter /ov e [0.017,0.018]. The tight 
constraints on log Dmix and /ov do not depend on the choice of 
the mass, metallicity, mixture and opacity table within the al¬ 
lowed seismic ranges. We also conclude that models with a step 
function core overshooting description are of inferior quality to 
explain the seismic data of this star. 

Our work represents the first application of seismic mod¬ 
elling of an SPB from which the need of global extra mixing is 
proven and quantified with high precision, based on the observed 
frequencies of unambiguously identified modes. Addressing the 
origin of the physical mechanisms that generate the desired value 
of L>mix, and their possible interaction(s), is yet to be done, keep¬ 
ing in mind that KIC 10526294 is an ultra-slow rotator. 

Forward seismic modelling is currently the best starting point 
to describe the global physical parameters, the interior struc¬ 
ture characteristics, and the detailed pulsational properties of real 
stars. To carry out an iterative asteroseismic inversion for the ro¬ 
tation profile or for the entire structure of the star, one has to start 
from a calibrated model that has the ability to fit all the detected 
and identified oscillation modes reasonably close to their mea¬ 
sured frequencies. In this study, we have achieved this for the 
slowly rotating Kepler target KIC 10526294. 

Despite the still relatively high values for the best 
model(s). Fig. 6 shows that the relative deviation between the 
modelled and observed frequencies is below 0.4%. Incorporating 
non-adiabatic corrections to the adiabatic frequencies, and mov¬ 


ing beyond the linearised oscillation paradigm (e.g. Van Hoolst 
1994) may further reduce the deviations between the observed 
and model frequencies. In addition, there is still room for im¬ 
proving the mode stability predictions for KIC 10526294 be¬ 
cause (a) depending on the adopted metallicity, several of the 
observed longest-period dipole modes are predicted to be stable, 
and (b) few excited dipole modes are not observed. 

Thanks to the unprecedented precision of the observations 
assembled by the nominal Kepler mission, we have now reached 
the stage at which we can employ asteroseismic methods to cali¬ 
brate the internal structure of massive stars with well-developed 
convective cores, and quantify the level of seismic modelling. 
This allows us to evaluate different choices of the input physics 
of B-star models coupled with an appropriate statistical frame¬ 
work (such as the x^q^ to properly account for the different of 
the models. Our modelling of KIC 10526294 proves that we can 
adequately explain the high-precision data collected from space 
within the framework of spherically symmetric stellar evolution 
models, and a linearised formalism of stellar oscillations below 
percent level. 
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Appendix A: Period spacing for iternative modeis 

The comparison between the period spacing of the first and sec¬ 
ond best models is presented in Fig. A.l. A similar plot for the 
second and third best models is shown in Fig. A.2. Both Model 8 
and Model 11 provide poorer fits to the period spacing on the 
short-period range of the series, than Model 4 (Fig. 4). 



Fig. A.l. Top panel: Period spacing for Model 4 (filled circles) and 
Model 11 (empty squares). See Table 3 for their parameters. Bottom 
panel: Absolute frequency deviations of the two models with respect 
to observations. See also Fig. 4 for a comparison. 



Fig. A.2. Similar to Fig. 4 but for Model 8 (filled circles) and Model 11 
(empty squares). 


Appendix B: Full overview of the Mixing, Fine 
and Metallicity Grids 

The Mixing, Fine and Metallicity grids introduced in Table 2 
have identical physical ingredients, and it is safe to merge 
their goodness-of-fit into a single snapshot. Fig. B.l shows 
log Ared ^ function of the dimensions of the three grids in addi¬ 
tion to the model radius and surface gravity. Note the finger-like 
structures in panels (c), (g), and (h) where the age, overshooting 
and mixing parameters are constrained, respectively. In panel (i), 
the position of all models is shown on the Kiel diagram (log Teff 
vs. log g), and that of the best model is marked with a (white) 
star. Our Model 4 is found marginally inside the 2cr uncertainty 
box from spectroscopy. 

Appendix C: Deliverables, inlists and opacity 
tables 

We adhere to the MESA code of conduct as stated in Paxton et al. 
(2011), and make our setting files and physical ingredients pub¬ 
licly available for download. This also ensures reproducibility of 
our results. The following items are available online: 

- The MESA v.5548 input inlists. 
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Fig. B.l. Behaviour of log;^^^^^ as a function of (a) effective temperature (b) surface gravity log g, (c) age Xc (d) radius, (e) initial mass, (f) 
metallicity Z, (g) exponential overshooting parameter /ov, and (h) extra diffusive mixing logD^ix. Panel (i) shows the position of the two best 
models (star: Model 4, pentagon: Model 11) in the Kiel diagram, where the box indicates the Icr and the dotted line the 3cr boundaries deduced 
from spectroscopy. The colour coding is based on log;^^^^^. 


The GYRE v.3.0 inlist, 

The OP and OPAL opacity tables adapted to the A05+Ne, 
A09, and NP12 mixtures. They have to be used with a 
standard MESA composition option initial_zfracs = 
5, 6 and 8, respectively. All tables are MESA compatible. 
When querying the OP and OPAL servers, we made a choice 
to redistribute the abundance residuals on all metals, based 
on their relative mass fraction instead of depositing them on 
the heaviest metals, which are Ee and Ni. This is a choice 


and may slightly influence the adiabatic and non-adiabatic 
results. 

- The internal structure of Model4, Model5, Models, and 
Model 11 in a GYRE-compatible format. 

Static links to download each of these products are available 
from https://fys.kuleuven.be/ster/Projects/ASAMBA. 
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